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INTRODUCTION 


This  annual  summary  report  describes  research  performed  from  1  July  1988  through 
31  July  1989  with  support  from  ONR.  During  the  last  ten  months  of  this  period 
(1  October  1988  through  31  July  1989),  the  support  came  from  ONR  grant  N00014- 
89-J-1003.  During  the  first  three  months  (1  July  1988  through  30  September  1988), 
the  support  was  provided  by  ONR  contract  N00014-85-K-0708. 

The  following  projects  are  discussed  in  this  report: 

I.  Scattering  of  Sound  by  Sound 

II.  Parametric  Receiving  Arrays 

III.  Nonlinear  Effects  in  Asymmetric  Sound  Beams 

IV.  Pulsed  Finite  Amplitude  Sound  Beams 

Contributions  to  these  projects  were  made  by  the  following  individuals: 

Senior  Personnel 

•  M.  F.  Hamilton,  principal  investigator  (projects  I— IV) 

•  J.  Naze  Tjptta,  visiting  scientist  (projects  I  and  II) 

•  S.  Tjptta,  visiting  scientist  (projects  I  and  II) 

Graduate  Students 

•  C.  M.  Darvennes,  Ph.D.  student  in  Mechanical  Engineering  (projects  I  and 

H) 

•  E.  E.  Kim,  M.S.  student  in  Mechanical  Engineering  (project  III) 

•  Y.-S.  Lee,  Ph.D.  student  in  Mechanical  Engineering  (project  IV) 

Naze  Tjptta  and  Tjptta,  who  spent  oast  year  at  Applied  Research  Laboratories 
at  the  University  of  Texas  at  Austin  '  ATL:UT)  while  on  leave  from  the  University 
of  Bergen,  Norway,  received  no  suppc.  *~om  the  ONR  contract  and  grant  covered 
by  this  report.  Their  support  on  the  above  projects  was  provided  by  ARL:UT, 
and  VISTA/STATOIL  of  Norway.  Kim  received  partial  support  from  the  National 
Science  Foundation  (2  months),  and  Lee  received  partial  support  from  the  Texas 
Advanced  Research  Program  (8  months). 

The  following  manuscripts  and  abstracts,  which  describe  work  supported  at 
least  in  part  by  ONR,  have  been  published  (or  submitted  for  publication)  since 
1  July  1988. 

Refereed  Publications 

•  M.  F.  Hamilton  and  J.  A.  TenCate,  “Finite  amplitude  sound  near  cutoff  in 
higher  order  modes  of  a  rectangular  duct,”  J.  Acoust.  Soc.  Am.  84,  pp.  327- 
334  (1988). 
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•  T.  S.  Hart  and  M.  F.  Hamilton,  “Nonlinear  effects  in  focused  sound  beams,” 
J.  Acoust.  Soc.  Am.  84,  pp.  1488-1496  (1988). 

•  C.  M.  Darvennes  and  M.  F.  Hamilton,  “Scattering  of  sound  by  sound  from 
two  Gaussian  beams,"  submitted  in  March  1989  for  publication  in  J.  Acoust . 
Soc.  Am. 

•  S.  J.  Lind  and  M.  F.  Hamilton,  “Noncollinear  interaction  of  a  tone  with 
noise,”  submitted  in  April  1989  for  publication  in  J.  Acoust.  Soc.  Am. 

Publications  in  Conference  Proceedings 

•  C.  M.  Darvennes,  M.  F.  Hamilton,  J.  Naze  Tjptta,  and  S.  Tjptta,  “Scattering 
of  sound  by  sound:  Effects  of  absorption,”  to  appear  in  Proceedings  of  the 
13th  International  Congress  on  Acoustics,  Belgrade,  Yugoslavia,  August  1989. 

Oral  Presentation  Abstracts 

•  Y.-S.  Lee  and  M.  F.  Hamilton,  “A  parametric  array  for  use  as  an  ultrasonic 
proximity  sensor  in  air,”  J.  Acoust.  Soc.  Am.  84,  S8  (1988). 

•  C.  M.  Darvennes,  M.  F.  Hamilton,  J.  Naze  Tj0tta,  and  S.  Tjptta,  “Effects  of 
absorption  on  the  scattering  of  sound  by  sound,”  J.  Acoust.  Soc.  Am.  85, 
S5  (1989). 

I.  SCATTERING  OF  SOUND  BY  SOUND 

The  effect  of  absorption  on  the  scattering  of  sound  by  sound  has  been  investigated 
theoretically  by  Darvennes,  Naze  Tjptta,  and  Tjptta.  The  contributions  from  Dar¬ 
vennes  are  closely  related  to  her  work  on  project  II.  Support  for  Naze  Tjptta  and 
Tj0tta  was  provided  by  ARL:UT  and  VISTA/STATOIL. 

A.  Background 

The  scattering  of  sound  by  sound  usually  refers  to  the  radiation  of  sum  or  difference 
frequency  sound  from  a  region  where  two  harmonic  sound  beams  interact  at  a 
nonzero  angle.  Background  on  investigations  into  the  scattering  of  sound  by  sound 
over  the  past  30  years  may  be  found  in  the  second1  and  third2  annual  summary 
reports  under  ONR  contract  N00014-85-K-0708. 

B.  Results 

The  results  obtained  during  the  past  year  are  summarized  by  the  following  abstract3 
of  an  oral  presentation  given  at  the  117th  Meeting  of  the  Acoustical  Society  of 
America  in  Syracuse,  New  York,  on  23  May  1989. 
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The  scattering  of  sound  by  sound  in  a  lossless  fluid  was  discussed  at 
an  earlier  meeting.4"’  Here,  the  effects  of  absorption  are  included.  The 
Khokhlov- Zabolotskaya-Kuznetsov  equation  is  used  to  derive  farfield 
asymptotic  results  for  the  sum  and  difference  frequency  sound  due  to 
the  noncollinear  interaction  of  real  sound  beams  radiated  from  displaced 
sources.  There  are  two  main  contributions  to  the  nonlinearly  generated 
sound  in  the  farfield:  the  continuously  pumped  sound  and  the  scattered 
sound.  Weak  absorption  affects  neither  the  locations  nor  the  relative 
amplitudes  of  the  pumped  and  scattered  difference  frequency  sound. 
Strong  absorption  attenuates  the  pumped  difference  frequency  sound 
faster  than  the  scattered  difference  frequency  sound.  The  scattered 
sum  frequency  sound  is  always  attenuated  faster  than  the  pumped  sum 
frequency  sound,  and  there  may  be  shifts  in  the  locations  of  the  maxima. 
Numerical  results  are  presented  for  the  case  of  Gaussian  primary  beams. 

The  numerical  results  that  follow  are  excerpted  from  a  paper6  to  be  published 
in  the  Proceedings  of  the  13th  International  Congress  on  Acoustics  to  be  held  in 
Belgrade,  Yugoslavia  during  August  1989.  We  consider  two  Gaussian  beams  that 
intersect  at  an  angle  of  10°.  The  sources  have  the  same  radius  a,  their  centers 
are  displaced  by  2a,  and  the  frequencies  of  the  two  beams  differ  by  a  factor  of 
five.  The  thermoviscous  absorption  of  the  fluid  is  characterized  in  terms  of  the 
absorption  length  L±  =  (cq  +a2-  a±)_1,  where  a  is  the  attenuation  coefficient  of 
the  primary,  sum,  or  difference  frequency  wave.  The  absorption  length  is  positive 
for  difference  frequency  generation  and  negative  for  sum  frequency  generation. 
Range  is  measured  in  terms  of  the  dimensionless  coordinate  Z  =  z/z2,  where  z2  is 
the  Rayleigh  distance  for  the  low  frequency  source.  The  dimensionless  size  of  the 
low  frequency  source  is  given  by  k2a  =  30,  where  k2  is  the  associated  wavenumber. 

In  Fig.  1  we  show  the  effects  of  moderate  absorption  on  sum  and  difference 
frequency  beam  patterns  at  ranges  from  Z  =  100  to  Z  =  500.  We  have  set 
\L±\  =  100*2,  and  therefore  the  dimensionless  range  Z  =  100  corresponds  to  one 
absorption  length.  In  the  difference  frequency  beam  patterns,  the  pumped  sound  is 
located  near  0  —  4.5°,  and  the  scattered  sound  is  located  near  0  =  7°.  Note  that  the 
shape  of  the  difference  frequency  beam  pattern  does  not  vary  with  range.  In  the 
sum  frequency  beam  patterns,  the  pumped  sound  is  again  located  near  6  =  4.5°, 
but  the  scattered  sound  is  now  located  near  0  =  3.5°.  Whereas  the  scattered  sum 
frequency  sound  dominates  the  pumped  sum  frequency  sound  at  Z  =  100,  the 
relative  levels  are  reversed  beyori  Z  =  400. 

In  Fig.  2  we  show  sum  and  difference  frequency  beam  patterns  at  Z  =  10  as 
the  absorption  length  is  varied  from  |L±|  =  oo  (no  absorption)  to  |L±j  =  z2  (strong 
absorption).  The  pumped  and  scattered  waves  are  located  near  the  same  angles  as 
in  Fig.  1.  Increasing  absorption  (decreasing  |L±|)  attenuates  the  pumped  difference 
frequency  sound  faster  than  the  scattered  difference  frequency  sound,  whereas  the 
reverse  is  true  for  the  sum  frequency  sound. 

We  conclude  that  the  relative  amplitudes  of  the  pumped  and  scattered  waves 
depend  not  only  on  the  nature  and  geometry  of  the  sources  (e.g.,  ka  of  each  source, 
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on-source  pressure  distributions,  source  separation,  and  beam  interaction  angle), 
but  also  on  the  dissipative  properties  of  the  medium  through  the  sign  of  L~±  — 
a(u,>!)  +  a(u> 2)  —  a(a>i  ±  u?).  The  function  a(a>)  describes  how  the  absorption 
coefficient  of  the  medium  varies  with  frequency.  At  distances  of  the  order  of  the 
absorption  length  |L±|  or  larger,  scattering  of  sound  by  sound  is  more  likely  to  be 
observed  at  combination  frequencies  for  which  L±  >  0. 

II.  PARAMETRIC  RECEIVING  ARRAYS 

The  performance  of  a  parametric  receiving  array  near  a  reflecting  surface  has  been 
investigated  theoretically  by  Darvennes  in  parallel  with  project  I.  Naze  Tjptta  (who 
is  a  member  of  Darvennes’  doctoral  committee)  and  Tj0tta  also  contributed  to  this 
work.  Research  into  this  problem  is  essentially  complete.  No  new  theoretical  work 
has  been  performed  since  the  publication  of  the  third  annual  summary  report2 
(July  1988)  for  ONR  contract  N00014-85-K-0708.  Darvennes  has  used  much  of  the 
last  year  to  report  her  work  in  a  dissertation  that  she  is  scheduled  to  defend  in 
August  1989. 
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A.  Background 

The  motivation  is  the  potential  use  of  a  parametric  receiving  array  for  measuring 
freefield  source  directivities  in  reverberant  environments.  Particular  attention  is 
devoted  to  sound  radiated  from  underwater  sources  near  the  air-water  interface. 
Background  for  this  work  may  be  found  in  the  second  annual  summary  report1  for 
ONR  contract  N00014-85-K-0708. 

B.  Results 

We  present  here  only  one  example  that  characterizes  the  results  from  this  project. 
A  comparison  is  made  between  two  different  methods  for  measuring  a  beam  pattern 
with  an  omnidirectional  hydrophone  in  the  farfield  of  an  underwater  sound  source. 
One  is  a  conventional  linear  measurement,  and  the  other  is  a  parametric  (nonlinear) 
measurement  that  makes  use  of  intermodulation  (difference  frequency)  components 
detected  by  the  hydrophone.  The  comparison  is  made  on  the  basis  of  which  method 
better  reproduces  the  beam  pattern  that  would  be  measured  in  the  absence  of  the 
reflecting  surface,  i.e.,  the  freefield  beam  pattern. 

As  an  example,  we  consider  an  axisymmetric  source  with  radius  10  m  that  radi¬ 
ates  in  water  at  250  Hz.  The  source  is  located  at  a  depth  of  50  m  below  the  surface, 
and  the  sound  is  radiated  in  a  direction  that  is  parallel  with  the  surface.  The  para¬ 
metric  array  is  formed  by  a  pump  with  radius  1  m  and  frequency  50  kHz.  The 
pump  is  located  on  the  low  frequency  source,  and  an  omnidirectional  hydrophone 
located  2.6  km  away  on  the  axis  of  the  pump  measures  not  only  the  primary  sound 
fields  at  250  Hz  and  50  kHz,  but  also  the  secondary  (difference  frequency)  sound 
field  at  49.75  kHz.  To  simplify  the  analysis  it  is  assumed  that  the  amplitude  shad¬ 
ing  of  both  the  source  and  the  pump  is  Gaussian,  and  that  the  water  is  infinitely 
deep. 

The  results  are  shown  in  Fig.  3.  The  dotted  line  (with  no  ripples)  is  the  freefield 
beam  pattern  of  the  250  Hz  source,  and  the  broken  line  is  the  beam  pattern  at 
250  Hz  that  would  be  measured  linearly  by  the  hydrophone  in  a  plane  that  is 
perpendicular  to  the  surface  of  the  water.  The  ripples  in  the  broken  line  result 
from  the  reflection  of  the  250  Hz  sound  from  the  surface  of  the  water.  The  solid 
line  is  determined  by  the  49.75  kHz  difference  frequency  sound  that  is  measured  by 
the  hydrophone.  Although  the  parametric  array  suppresses  the  ripples  caused  by 
the  reflecting  surface,  it  tends  to  overestimate  the  beamwidth  of  the  low  frequency 
source.  Moreover,  there  is  an  optimal  range  for  the  separation  between  the  pump 
and  the  hydrophone.  As  the  separation  is  increased,  the  signal  from  the  parametric 
array  is  increasingly  affected  by  the  multipath  components,  and  as  the  separation  is 
decreased,  the  beamwidth  of  the  low  frequency  source  is  increasingly  overestimated. 
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Figure  3 


III.  NONLINEAR  EFFECTS  IN  ASYMMETRIC  SOUND 
BEAMS 

This  is  a  new  project  that  was  begun  in  July  1988  by  Kiin.  Kim  received  two 
months  of  support  from  the  National  Science  Foundation  (1  July  1988  through 
31  August  1988)  during  the  period  covered  by  this  report. 

A.  Background 

The  effects  of  asymmetry  on  nonlinear  distortion  are  relevant  to  bottom  and  sub¬ 
bottom  profiling  with  a  parametric  array,  high  intensity  sound  beams  in  shallow 
water,  and  caustic  formation  in  inhomogeneous  media.  In  general,  asymmetries 
are  introduced  whenever  a  sound  beam  encounters  sin  interface  from  a  direction 
that  deviates  from  the  normal.  For  example,  Muir  et  al.7  have  shown  that  near 
critical  incidence,  a  parametric  array  attains  greater  bottom  penetration  through 
a  water-sediment  interface  than  does  a  comparable  beam  generated  by  a  linear 
source.  Theoretical  explanations  of  this  phenomenon  have  been  offered  by  Jarzyn- 
ski  and  Flax8  and  by  Berktay  et  al.,9  although  diffraction  effects  are  excluded  in 
both  analyses.10’11,12  Only  recently  have  appropriate  theoretical  models  .*or  diffract- 
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ing  finite  amplitude  sound  beams  at  interfaces  and  in  inhomogeneous  media  been 
developed  by  the  Tjottas  and  their  doctoral  students  (e.g.,  Ref.  13). 

Our  research  is  limited  to  two  dimensional  sound  fields  because  of  the  rel¬ 
ative  ease  with  which  numerical  computations  can  be  made  at  high  intensities 
(in  comparison  with  the  same  calculations  for  three  dimensional  nonaxisyinmetrie 
fields).  Investigations  based  on  two  dimensional  fields  are  nevertheless  appropriate 
for  waveguides,  and  also  for  comparison  with  recent  analytical  work  on  caustics 
by  Marston.141"'  McDonald  and  Kuperman16  have  recently  produced  numerical  re¬ 
sults  for  the  water-borne  propagation  of  a  two  dimensional  finite  amplitude  pulse 
through  a  caustic. 

B.  Results 

The  work  is  based  on  a  computer  program  for  axisymmetric  sound  beams  that  was 
developed  originally  by  Aanonsen  et  al. 17,18  References  19  and  20  discuss  subse¬ 
quent  modifications  of  the  program.  The  program  has  been  modified  by  Kim  for 
application  to  two  dimensional  asymmetric  sound  beams. 

Shown  in  Fig.  4  are  preliminary  results  for  the  beam  patterns  in  a  high  inten¬ 
sity  sound  beam  radiated  by  an  infinite  strip  source.  The  solid  lines  are  for  the 
fundamental  component,  and  the  successively  lower  broken  lines  correspond  to  the 
nonlinearly  generated  higher  harmonic  components  (second  through  fourth).  The 
source  excitation  is  time  harmonic,  and  the  pressure  amplitude  decreases  linearly 
from  po  at  one  edge  to  (1  —  e)po  it  the  other,  as  shown  schematically  underneath 
each  set  of  beam  patterns.  Source  asymmetry  is  thus  characterized  by  the  dimen 
sionless  parameter  e.  With  e  =  0  the  source  has  a  uniform  amplitude  distribution, 
and  the  resulting  beam  patterns  resemble  those  for  the  three  dimensional  case  of 
a  baffled  piston.  The  effects  of  asymmetry  at  e  =  0.5  are  most  prominently  mani¬ 
fested  by  the  nonlinear  nearfield  phenomena  known  as  fingers21  (i.e.,  the  additional 
sidelobes  in  the  second,  third,  and  fourth  harmonic  beam  patterns).  The  asym¬ 
metry  of  fhe  fingers  increases  with  the  number  of  the  harmonic.  At  e  =  1,  the 
fingers  can  no  longer  be  observed,  and  the  overall  shapes  of  the  beam  patters  are 
approximately  symmetrical. 

IV.  PULSED  FINITE  AMPLITUDE  SOUND  BEAMS 

This  is  a  new  project  that  was  begun  by  Lee  in  September  1988.  Lee  received  eight 
months  of  support  from  the  Texas  Advanced  Research  Program  ( 1  September  1988 
through  30  April  1989)  during  the  period  covered  by  this  report. 

A.  Background 

Although  a  broad  theoretical  base  exists  for  the  linear  analysis  of  transient  radi¬ 
ation  in  directive  sound  beams,22,23  the  investigation  of  transient  effects  in  finite 
amplitude  sound  beams  has  progressed  very  slowly  over  the  years.  The  first  the¬ 
oretical  model  of  transient  radiation  from  a  directive  finite  amplitude  source  was 


7 


developed  by  Berkt avJ  1  soon  after  Westervelt  s  discovery  of  the  pararnet ric  array.2"' 
Berktay's  model  applies  to  a  finite  amplitude  sigr  \1  whose  source  waveform  is 
a  slowly  modulated,  high  frequency  carrier  wave.  Two  simplifying  assumptions 
are  inherent  in  his  analysis:  (1)  the  primary  beam  behaves  as  a  nondiffracting, 
collimated  plane  wave  whose  finite  amplitude  distortion  is  terminated  within  the 
nearfield  by  absorption,  and  (2)  the  results  are  valid  only  in  the  farfield.  However, 
Berktay’s  prediction  that  the  farfield  time  waveform  should  be  described  on  axis 
by  the  second  time  derivative  of  the  square  of  the  source  envelope  function  seemed 
to  be  verified  by  the  classic  experimental  work  of  Moffett  and  coworkers.26,27,28 

Subsequent  extensions  of  Berktay’s  original  analysis  have  been  few  and  far  be¬ 
tween.  Singhal  and  Zornig2!>  considered  the  inverse  problem,  the  calculation  of  the 
source  waveform  required  to  obtain  a  desired  farfield  response.  Trivett  and  Rogers30 
investigated  theoretically  the  second  order  sound  that  is  generated  by  the  nonlinear 
interaction  of  a  pulsed  beam  which  intersects  the  path  of  a  monofreq  lency  plane 
wave  at  an  arbitrary  angle.  In  a  later  paper31  the  authors  extended  their  analysis  to 
include  the  case  where  the  monofrequency  wave  is  a  collimated  plane  wave  beam. 
The  collinear  interaction  of  nondiffracting,  collimated  plane  w'ave  pulses  was  also 
investigated  by  Pace  and  Ceen32  and,  more  recently,  by  Stepanishen  and  Koenigs.3,3 
Common  to  each  of  these  analyses,  however,  are  the  same  restrictions  which  apply 
to  Berktay’s  work,  namely,  that  the  primary  waves  experience  no  diffraction,  that 
the  results  apply  only  to  the  farfield  radiation,  and  that  the  beam  exhibits  only 
moderately  nonlinear  effects. 

The  first  theoretical  analyses  to  take  into  account  the  diffraction  of  pulsed  fi¬ 
nite  amplitude  sound  beams  were  performed  in  the  Soviet  Union.  Zabolotskaya  and 
Khokhlov,34  in  their  paper  where  they  derive  the  lossless  nonlinear  parabolic  wave 
equation  which  bears  their  names,  present  a  theoretical  analysis  of  simple  unipolar 
pulses  in  diffracting  sound  beams.  Zhileikin  and  Rudenko,35  using  a  numerical  so¬ 
lution  of  the  Khokhlov-Zabolotskaya  equation,36  investigated  unipolar  and  bipolar 
Gaussian  pulses  in  finite  amplitude  Gaussian  beams.  However,  the  more  practical 
case  of  tone  bursts  radiated  from  piston-like  sources  was  not  considered. 

Frpysa,  Naze  Tjotta,  and  Tjptta37,38  recently  developed  a  quasilinear  theory 
that  takes  into  account  the  combined  effects  of  diffraction,  absorption,  and  moder¬ 
ate  nonlinearity  in  pulsed  sound  beams.  Their  work  shows  that  the  classic  result 
by  Berktay,24  which  seemed  to  be  confirmed  experimentally,26,27’28  becomes  increas¬ 
ingly  inaccurate  as  the  number  of  cycles  in  the  pulse  is  reduced  (i.e.,  for  short  tone 
bursts).  The  apparent  experimental  confirmation  of  Berktay’s  theory  was  due  to 
the  fact  that  relatively  long  tone  bursts  were  used.  Results  from  the  quasilinear 
theory  will  provide  useful  comparisons  with  the  results  obtained  from  the  algorithm 
being  developed  by  Lee  for  strong  finite  amplitude  pulses. 

To  date,  there  are  no  analytical  solutions  that  accurately  describe  the  com¬ 
bined  effects  of  diffraction,  absorption,  and  strong  nonlinearity  in  sound  beams. 
Instead,  numerical  solutions  are  used.  The  most  common  model  equation  for  finite 
amplitude  sound  beams  is  the  Khokhlov-Zabolotskaya-Kuznetsov  (KZK)  nonlin¬ 
ear  parabolic  equation.  An  algorithm  developed  by  Aanonsen17  that  solves  the 
KZK  equation  in  the  frequency  domain  has  been  used  with  great  success  for 
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monofrequency16,19,39-40  and  bifrequency41'42  sources.  However,  the  number  of  har¬ 
monics  that  must  be  retained  to  describe  the  distortion  in  pulsed  sound  beams 
may  require  prohibitive  amounts  of  computer  time.  Bj0rn0  and  Neighbors43  have 
recently  developed  a  scheme  by  which  the  Aanonsen  algorithm  may  be  used  for  ccr 
tain  pulsed  sound  beams.  To  reduce  the  number  of  harmonics  required  to  describe 
a  single  pulse,  they  consider  an  infinite  series  of  pulses  that  are  spaced  sufficiently 
far  apart  in  time  that  adjacent  pulses  do  not  interact  with  each  other.  Neverthe¬ 
less,  pulses  with  short  rise  times  are  inherently  difficult  to  describe  numerically  in 
the  frequency  domain,  particularly  when  nonlinear  effects  are  relatively  strong. 

An  alternative  approach  is  provided  by  a  time  domain  algorithm  that  was  devel¬ 
oped  by  McDonald  and  Kuperman. 16,44  Their  program  solves  a  nonlinear  parabolic 
wave  equation  that  is  similar  to  a  KZK  equation  without  the  absorption  term. 
However,  their  method  for  integrating  the  diffraction  term  can  become  unstable  in 
the  highly  oscillatory  nearfield  of  a  baffled  piston. 

B.  Results 

Lee  has  developed  a  numerical  algorithm  that  solves  the  KZK  equation  in  the 
time  domain.  Very  near  the  source,  where  the  solution  can  oscillate  very  rapidly, 
the  diffraction  term  is  integrated  separately  from  the  dissipation  and  nonlinearity 
terms.  The  spatial  integration  of  the  diffraction  term  is  accomplished  with  an 
implicit  backward  finite  difference  method  (as  in  the  Aanonsen  algorithm),  and  a 
trapezoidal  rule  is  used  to  integrate  over  time.  The  nonlinearity  and  dissipation 
terms  are  integrated  simultaneously  with  a  predictor-corrector  method,  as  has  been 
done  for  the  Burgers  equation.45  Away  from  the  source,  a  transformed  version  of 
the  KZK  equation  is  used,  and  all  terms  are  integrated  simultaneously  with  an 
alternating  direction  implicit  method. 

Lee  is  still  in  the  process  of  checking  his  numerical  results  against  exact  ana¬ 
lytical  solutions  for  various  limiting  cases.  Results  from  numerical  integration  of 
the  dissipation  and  nonlinearity  terms  (i.e.,  the  Burgers  equation)  are  in  excellent 
agreement  with  the  Hopf-Cole  analytical  solution.  Shown  in  Fig.  5  are  results  for 
the  diffraction  term  alone,  which  presented  the  most  difficulty  for  numerical  inte¬ 
gration.  The  case  considered  is  the  on-axis  field  of  a  baffled  piston  that  radiates 
a  single  cycle  of  a  sine  wave,  as  shown  in  the  top  two  figures.  The  left  column 
shows  the  numerical  solution  and  the  right  column  shows  the  analytical  solution. 
The  dimensionless  coordinate  a  measures  range  in  terms  of  the  Rayleigh  distance 
for  an  infinite  sine  wave.  Beyond  a  «  0.3,  the  numerical  solution  is  in  excellent 
agreement  with  the  analytical  solution,  apart  from  a  small  numerical  disturbance 
that  appears  at  the  trailing  edge  of  the  pulse.  Solutions  for  the  highly  oscillatory 
field  in  the  neighborhood  of  baffled  piston  sources  are  also  difficult  to  obtain  when 
•  e  KLK  equation  is  integrated  numerically  in  the  frequency  domain.18 
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